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ABSTRACT 

When a gluon or a quark is sent through the hot QCD plasma it can be absorbed into the 
ambient heat bath and so can acquire an effective lifetime. At high temperatures and for 
weak couplings the inverse lifetime, or damping rate, for energetic quarks and transverse 
gluons, (those whose momenta satisfy |p| ^ gT) is given by 7(p) = c (7^ log T + 
0{g'^T). We show that very simple arguments suffice both to fix the numerical coefficient, c, 
in this expression and to show that the 0{g'^T) contribution is incalculable in perturbation 
theory without further assumptions. For QCD with Nc colours we find (expressed in terms 
of the casimir invariants Ca = Nc and Cf = {N^ — l)/{2Nc)): Cg = for gluons and 
Cq = +^ for quarks. These numbers agree with the more detailed calculations of Pisarski 
et.al. but disagree with those of Lebedev and Smilga. The simplicity of the calculation 
also permits a direct verification of the gauge-invariance and physical sign of the result. 
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The behaviour of nonabehan gauge theories at finite temperature is of theoretical 
interest due to the surprisingly rich structure they exhibit even within the perturbative 
regime of weak coupling and high temperatures. They may also have phenomenological 
applications to the interpretation of the collisions of heavy nuclei at high energies. 

In recent years much attention has been directed towards understanding the properties 
that gluons and quarks acquire as they propagate through a quark-gluon plasma. -"^ This 
has been at least partially due to the sometimes contradictory and confusing results of the 
earliest one-loop calculations. In these calculations the gluon damping constant, 7g(p), 
was found in the static (|p| — > 0) limit to depend on the gauge on which it was computed 
and, for some gauges, to be negative. If taken seriously such a negative damping constant 
would indicate an instability of the thermal state towards gluon emission. The source of 
the confusion with these calculations^ is that they neglect higher-loop contributions that 
are of the same order in the gauge coupling constant, as are the terms that are kept. 
This is because in perturbation theory (in three space dimensions) at finite temperatures 
successive terms in the loop expansion need not be suppressed relative to fewer loops by 
additional factors of the gauge coupling. 

The failure of the loop expansion arises at finite temperatures because of the occurence 
of severe infrared divergences. These infrared divergences are more troublesome than they 
are at zero temperature due to the singular behaviour of the Bose-Einstein distribution 

function at low energies: n{E) = {e^l^ — k, ^ -| . This behaviour causes quantities 

to blow up like a power of the infrared cutoff rather than simply logarithmically as they do 
at zero temperature.^ Indeed, a simple power-counting argument shows that for QCD at 
temperature T a generic £-loop graph can contribute an amount proportional to {g'^T/XY 
relative to the tree-level result. A in this expression is an infrared cutoff. Higher loops 
clearly need not be suppressed once A is as low as g'^T. 

In fact the loop expansion already fails even for an infrared cutoff as large as A ft! gT 
since a subclass of diagrams can be even more infrared singular than is indicated by the 
generic power-counting argument. The dangerous graphs are those such as the 'ring' graphs 
within which multiple self-energy insertions are made along a single internal line.^ In 
reference [5] Braaten and Pisarski argue that these last contributions may be resummed by 
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dressing all 'soft' lines — those carrying momenta less than or of order gT — by the calculable 
contributions of 'hard thermal loops'. In ref. [2] they then compute the implications 
for the particular case of the damping rate for static gluons, arguing that the failure 
of perturbation theory is in this case completely cured by such a resummation. This 
last conclusion is, however, disputed by Lebedev and Smilga,^ who argue that other 
contributions beyond the hard thermal loops can contribute equally large effects to the 
static gluon damping constant. 

The purpose of this letter is to point out that for at least some physical quantities 
the dominant part of the result for small g can be identified quite simply without making 
use of the complete resummation formalism. We take by way of illustration the damping 
constant, 7(p), for transverse gluons (and for quarks) but evaluated for momenta |p| » gT 
rather than in the static limit. For these momenta the dominant contribution for small 
coupling to 7(p) is momentum independent and has the form c g^ log j T + 6 g^T + ■ ■ ■. 
We show that the coefficient c can be computed by a simple, analytic calculation that 
only involves wavelengths that are within the perturbative regime: A ^ g^T. The same is 
not true for the coefficient h which can receive contributions for A ^ g^T. We can in this 
instance therefore explicitly verify that 7(p) is indeed gauge- independent and positive. 

Our results, given in and immediately following eq. (12), may be compared to more 
detailed calculations. For fermions they agree with those of ref. [7] but disagree by a 
factor of 3 with the real-time calculations of ref. [6] . For transverse gluons we agree with 
the result of an as-yet-unpublished version of a full resummation calculation,® but do not 
agree with the real-time estimate of ref. [6] . 

We now turn to a description of our calculation. We focus first on the purely gluonic 
theory since this is the case for which there is the most simplification over a full calculation. 

The dispersion relation relating the energy and momentum for a relativistic transverse 
gluon traversing the plasma may be determined by the position, -E(p) = — ■17 (p), in 
the complex energy plane of the zero of the transverse part of the inverse gluon propaga- 
tor. Using the usual one-parameter family of covariant gauges for which the bare inverse 
Feynman propagator is: 

iG^L),. ip) = P'ri,. + (r' - 1) P,P.. (1) 
3 



and the full thermal inverse propagator is: 



(GTun),. ip) = {G^D,. iP) + ^,AE, P), (2) 
the dispersion relation for transverse gluons becomes: 

(3) 



n* - 



Here i and j are to be summed over the three spatial directions, i,j = 1,2,3, in the 
plasma rest frame. Our task is to compute the leading contributions to Il^i,{E, p) and to, 
in particular, identify the origin of the log g terms. 

In order to sort out the size of the various contributions to Il^i, it is useful to divide 
the loop integrations according to whether or not they involve only energies and momenta 
that are greater than some infrared cutoff, A. That is: 

U^^iE, p) = UlfiE, p. A) + Ull^'^iE, p. A), (4) 

in which all loop momenta in 11^^'^^ {E, p. A) are cut off in the infrared at A. Provided that 
A is chosen sufficiently high H^^^^ is calculable within the loop expansion. 11^°^ is not 
computable in this way and must be obtained by other means. An important issue to be 
addressed is how much any desired quantity depends on the largely undetermined n^*^^*. 

If A is chosen much greater than gT then the lowest-order contribution to H^^^ arises 
at 0{g'^) due to one-loop graphs. There are four such one-loop graphs: one tadpole and 
one vacuum-polarization graph having either an internal gluon, quark or ghost loop. For 
example, using the Matsubara imaginary-time technique,^ the gluon vacuum-polarization 
graph contributes an amount: 

2/^ rp ^ f 73 

n;i-d(£;,p,A) = ^-^ ^ / ^v;Ap(p,-P-?,?)G^np + ?) 

X G'^'^(g) Ka«(-P,-g,P + g)- 

Here the gluon propagator is the inverse of eq. (1) and the three-point gluon vertex is given 

by the zero-temperature expression: V^\p{p, q, k) = {p-q)pr]^x\ + («? - /c)^??pA + {k-p)\r]pp. 
The quadratic invariant in the adjoint representation is Ca = Nc (= 3) for the gauge group 



SU{Nc). We work in Euclidean signature and the time component of every four- vector 
is an integer times 27rT. The summation is over the integer corresponding to the loop 
momentum, q. The lower limit on the momentum integration is meant as a reminder of 
the infrared cutoff. 

If the infrared cutoff should instead be chosen to satisfy g'^T <^ A <^ gT then higher- 
loop graphs can contribute to the same order in g as do the one-loop graphs in which 
the loop momenta are of order gT. After resumming these higher-loop contributions 
the leading result for nj^^'^'^ is given by the same one-loop graphs as before, with the 
proviso^ that each propagator (or vertex) is to be replaced with an 'effective' resummed 
propagator (or vertex): G^j^ — > G*^, and F^Ap — ^^Xp- Each of these resummed quantities 
agrees (up to higher powers of g) with its bare counterpart unless all of the momenta 
entering the line (or vertex) in question are themselves 0{gT). For such soft momenta, 
however, the resummed items are relatively complicated functions of the ratios of the 
energies and momenta that flow through that part of the graph. The main technical 
difficulty with working with the effective propagators and vertices lies in manipulating this 
more complicated form. 

Some of this extra complication simplifies if the external momentum is itself chosen to 
be much larger than gT, as is the case for when the dispersion relation, E{p), is evaluated 
for hard momenta such as |p| T. The reason for this simplification is that since the 
resummed and bare quantities agree if any of the relevant momenta are hard, |p| ^ gT, it 
is only necessary to work with the resummed versions if all of the momenta passing through 
a particular line (or vertex) are soft. For hard external momenta many of the internal lines 
and vertices must also carry hard momenta and so may be represented by the usual bare 
Feynman rules. For example, in the vacuum-polarization graphs both of the vertices and 
one propagator must necessarily carry hard momenta if the external momentum is itself 
large in comparison to gT. This leaves at most a single internal propagator to be dressed. 

Now comes the main point: the dominant part of H^j^ for small coupling — i.e. the 
g'^ logg terms — can be determined with no knowledge of H^'J,^*, and using virtually none of 
the complications of the resummation formalism in 11^^'^. 

The principal observation is that all of the terms in n^i^ that are proportional to 
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^^log^ are completely determined by the infrared divergent part of nj^^'^'^. To see how 
this works consider the lowest order contribution to the dispersion relation, -E(p), for 
transverse gluons with momenta |p| ~ T ^ gT. This is determined from eq. (3) given the 
vacuum polarization, Il^^{E,p) evaluated at hard, on-shell momenta: E — \p\ T. As 
is established in more detail below, the contribution of H^^^ to the right-hand-side of eq. 
(3) diverges logarithmically with A for g'^T <^ A <S gT: 



F'^^'^^(r,|p|,A) = 



1 



hardy _^ ^ »i 

/ i 



= 9 



(n 

^log 



A 



A^hard 



+ B + 



(6) 



+ Oig'). 



In this expression A and B are purely functions of |p| and T and //hard is a calculable energy 
scale of the high-frequency part of the theory which turns out below to be Hhard ~ gT. 

In order to extract information about F(T, |p|) = F^''"^{T, |p|, A) + F"°f^(T, |p|,A) 
from eq. (6) it is necessary to say something about the behaviour of F^°^^{T, |p|, A). The 
only property of F^°^^{T, |p|, A) that is required is that its A— dependence must cancel that 
of F^^'^(T,\p\,X): 



A — :rT = — A 



dX 
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dX 
-A + 



(7) 



+ 0{g'). 



This determines F^°^^ to have the form: 



psoft 
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log(^) +C + o(^ 



V A J 



+ 0{g'), 



(8) 



in which A is the same function as in eq. (6). The constant //soft that appears within the 
logarithm in this equation is the constant of integration that arises in passing from eq. (7) 
to eq. (8). It has dimensions of mass and is chosen to be of order g'^T since this is the 
largest mass scale present in the soft part of the problem. 
Adding the results of eqs. (6) and (8) therefore gives: 



F = -Ag^ log +g\B + C)+0 + 0{g') 



-Ag'log[^)+0{g'), 



(9) 
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which determines the coefficient of the log g term completely in terms of the calculable 
coefficient A. 

The next point is that since A is determined by the infrared divergent part of F^^'^'^^ it 
is insensitive to most of the complications of the resummation. To illustrate the simplicity 
with which A may be determined we now outline its calculation for the damping constant. 

Inspection of the Feynman rules shows that the only potentially infrared-divergent 
part of the leading contribution to F^^^'^ comes purely from the vacuum-polarization graph 
for nj^^'^'^ in which it is a gluon which circulates around the loop. Furthermore, even 
for this graph an infrared divergence can arise only from the term for which the integer 
n = q^/ (27rT) for the soft gluon line in eq. (5) vanishes, and even then only if the external 
four-momentum is on shell: E = |p|. Since only n = contributes, it is sufficient to know 
the form for the resummed propagator at zero frequency, where it reduces to:^ 

[(GT'],. ip) - {G^D,. iP) + m'P,.. (10) 

Here P^y is the projection matrix onto the rest frame of the plasma, and so is given in 
this frame by the matrix diag(0, 1, 1, 1). m denotes the lowest order gluon mass, or plasma 
frequency, which is given in terms of the number of quarks, n^, by w? = ^{gT)"^ (Ca + ^) . 

Substituting this into eq. (5), evaluating on the lowest-order mass shell, E = |p| fa T, 
and recognizing that at most one internal line can be soft at a time (and so need be dressed) 
then gives an infrared-divergent contribution: 



(^'|Pl)lciiv = 2^], q{q-+W "^ 

_,^g^C^T\p\^fX 



q-2\ 


IpI 




g + 2|p| 





2tt \m 



(11) 



from which we read A = (zCoT|p|)/(27r) and //goft = rn. Notice that the real part of F is 
infrared finite so only the lifetime acquires a g^ log g contribution. Using this result in the 
mass-shell condition gives our main result: 

ImF 



7<?(P) 



2|p| 



^^«^log(l)+0(^^). 



47r \g 
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(12) 



This entire argument when repeated for the quark self-energy similarly gives a log g 
contribution to the branch of the fermion spectrum that survives at large momenta. A 
simple calculation gives its coefficient as 7g(p) = ^ log j + 0{g'^). Cf denotes the 
quadratic invariant in the fundamental representation: Cf — {N^ — l)/{2Nc) (= |). 

There are several features of this calculation that bear emphasis: (i) First, as is re- 
quired for a good approximation to a physical quantity, 7 is independent of the gauge pa- 
rameter ^. All terms that depend on the gauge-parameter contribute only an infrared-finite 
result, (a) The sign of 7 is positive, indicating stability, (iii) Notice that a determination 
of the subleading 0{g'^) contributions would require knowledge of both of the coefficients 
B of eq. (6) and C of eq. (8). Although B is calculable using the complete resummation 
formulation, C is not and can only at present be determined by making some assumptions 
concerning the behaviour of the plasma in the low-frequency regime A ^ g'^T. It follows 
that the coefficient B need not a priori by itself be gauge-independent (or positive), (iv) 
Also, since it is a logarithmic infrared-divergence that is responsible for the logarithmic de- 
pendence on g, its coefficient is insensitive to the details of how the cutoff is implemented. 
Finally, (v) since the infrared-divergent term in F^ard explicitly proportional to {cf. 
eq. (11)) it only receives contributions from the two-loop and higher graphs that serve 
to dress the soft propagator in the gluon vacuum-polarization graph. It also follows that 
the imaginary part arises only from the self-energy of internal lines which carry soft loop 
momenta |q| < gT since the mass m may be taken to be zero for larger momenta. This 
agrees with what is expected physically from unitarity given the constraints of energy and 
momentum conservation in the plasma. 

We conclude that for some quantities in which infrared divergences in the perturbative 
expansion introduce a logarithmic dependence on the gauge coupling, g, it is possible to 
very simply identify the dominant contributions. This simplicity allows a check on more 
complete and more involved calculations. 
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